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The approximation of parameterized curves by segments of parabolas that pass 
through the endpoints of each curve segment arises naturally in all quadratic iso- 
parametric transformations. While not as popular as cubics in curve design problems, 
the use of parabolas allows the introduction of a geometric measure of the discrepancy 
between given and approximating curves. The free parameters of the parabola may be 
used to optimize the fit, and constraints that prevent overspill and curve degeneracy 
are introduced. This leads to a constrained optimization problem in two variables 
that can be solved quickly and reliably by a simple method that takes advantage of the 
special structure of the problem. 


FORMULATION OF THE PROBLEM 


We assume that tessellation has been done in such a way that every curve segment 
lies to one side of the straight line segment that connects consecutive nodes, and 
that the approximating curve segment results from a quadratic isoparametric transfor- 
mation and is therefore parabolic [1] and has the nodes as endpoints. The free 
parameters of the transformation are the coordinates of that point (a, 3) on the 
approximating parabola where the slope is identical to that of the node-connecting 
straight-line segment . We impose constraints that prevent image overspill (equiva- 
lently, vanishing of the transformation Jacobian over an associated triangular master 
element) and ensure limit stability as the curve degenerates to a straight line. We 
introduce coordinate axes in such a way that the nodes are at (-1,0) and (1,0), with 
the curve segment above the horizontal x-axis. The approximating parabolic segments 
can then be expressed in parameterized form as 

x = t + a(l -t 2 ), y = 8(1 -t 2 ), -1 < t < 1 (1) 


where the parameters a and 3 completely determine a particular member of the approxi- 
mating set (1). The axis of the parabola has slope — . The constraints associated 
with overspill have the form 

| - a — L 3 > o, |+a — — B > 0 (2) 

y l y 2 


where and y 2 are constants that are related to the transformation [2], We ensure 
limit stability by restricting the axial slope [1], namely 



(3) 


where m is a positive stability safety factor. The constraints (2) and (3) define the 
feasible domain in which the point (a, 3) may lie. 

The parameters a and 3 may be used to optimize the approximation of a given curve 
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x = A(s) , y = B(s), - 1 < s < 1 


(4) 


with 


A ( 1 ) = - A(-l) = 1, B ( 1 ) = B(-l) = 0, B (s) > 0 for - 1 < s < 1 

once we have defined a discrepancy measure for (1) and (4). When the approximating 
curves are parabolas, a uniquely defined distance between any point on the given 
curve and a parabola can be measured in the direction of the axis of the parabola. 

This results, in effect, in a re-parameterization of the parabola with respect to the 
given curve parameter s. Integration of the squares of the distances for the relevant 
parameter interval - 1 < s < 1 yields the objective function 

2 1 

F(a,3) = (\ + D / (3(1 - t 2 (s)) - B(s)) 2 ds (5) 

3 -1 

where 

t(s) = A(s) - | B(s) 

This paper is concerned with the optimization of (5) subject to the constraints of (2) 
and (3). 


THE OPTIMIZATION PROCEDURE 


a 


When the inverse g- of the axial slope is considered as a new variable y, (2), 

(3) and (5) can be expressed in terms of 3 and y, and the equivalent constrained 
optimization problem can be formulated, namely, minimize 


g(3,y) = (Y 2 + D(3 2 c 2 (y) + 3C 1 (y) + c 0 ) 


( 6 ) 


where 


1 9 9 1 9 

c (y) = / (1 - t (s)) ds , C (Y) = - 2 7 B(s) (1 - t: (s))ds 


-1 


-1 


C n = / B 2 (s)ds , t(s) = A(s) - YB(s) 

-1 


'0 


(7) 


subject to 


3 > 0, - | < Y s ^ , 5 - 3(Y + 7 ~) > 0, i + 3(Y -TT-) > 0 


1_ 

m 


( 8 ) 


The main advantage of the formulation (6) to (8) is that we can now use a separated 
optimization procedure [3] to find the best fit. We define the relevant objective 
function 


g(y) 


min 


0 * 3 * 3 . 


max 


g(3,y) 


(9) 
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where 3 is a function of y, and 3 is restricted to the feasible domain by the 

max ^ 11 
upper limit. The function G(y) is optimized with respect to y, for - — < y < — , to 

find the optimal value y. The optimal 3 is obtained by introducing y in (9) and 
optimizing with respect to 3- Justification of this procedure is presented in [4], 
and its implementation is described in [5]. It is reliable, efficient and eminently 
suitable for use on smaller computers and in the absence of multi-purpose optimization 
software. 


DISCUSSION OF APPLICATIONS 


For applications in the field of computer-aided design, the given curves (4) are 
often cubic polynomials, and the coefficients (7) may be calculated in closed form in 
terms of the polynomial coefficients by using a symbolic machine language so that 
families of curves can be approximated with no further integration. For general 
curves, numerical quadrature may be used, as in the implementation [5] where the Romberg 
quadrature is applied. The coefficient functions C^ (y) and C^Cy) are expanded as 

polynomials in y, so that for given A(s) and B(s) the integrations need only be done 
once. 


The method was used to find optimal constrained parabolic approximation to a wide 
variety of given curves. Some examples from [6] were included in the numerical tests. 
A comprehensive discussion of the experimental results is contained in [4]. The 
method yielded satisfactory approximations to the given curves for all the examples 
considered. 
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